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PREFACE 

This pamphlet, one of the series which will c >mprise the Ordnance Engineering 
Design Handbook, deals with the subject of Exterior Ballistics, covering Tra- 
jectories, Differential Effects and Data for Projectiles. The pamphlet pertains 
to that phase of projectile life starting at the gun muzzle (or, for a rocket, at 
burnout ) and ending at the point of burst or impact. The subject is covered in a 
manner intended primarily to assist the ammunition designer in recognizing and 
dealing with the design parameters of projectile weight, shape, muzzle velocity, yaw, 
drag and stability, in designing for optimum results in range, time of flight and 
accuracy. It is hoped that it will also prove helpful to other personnel of the Ord- 
nance Corps and of contractors to the Ordnance Corps. Because of its complexity, 
it is possible to cover the subject only briefly in this volume, but an attempt has been 
made to incorporate references which will permit the designer or the student to 
explore in more detail any phase of the subject. 

The ballistics of bombs and of intermediate and long range missiles are not 
covered in this pamphlet. Neither are the problems associated with Interior 
Ballistics or of Terminal Ballistics. These subjects will be covered by other 
pamphlets. 
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SYMBOLS 



A — Azimuth, measured from the direction of mo- 
tion of the airplane (in aircraft fire) 

A — Altitude function (Eq. 63) 

A — Axial moment of inertia 
B — Drag coefficient in lb/in s -ft (Eq, 19) 

C — Ballistic coefficient (Eq. 22) 

D — Drag 

D — Angular deflection 
F — Temperature in degrees Fahrenheit 
G — Drag function (Eq. 20) 

H — Ratio of air density at altitude y to that at 
the surface (Eq. 48) 

H — Height of target 
I — Inclination function (Eq. 62) 

K — Coefficient ( See Eq. 16 and 34) 

L — Crosswind force (lift) 

M — Mach number (Eq. 18) 

M — Pseudo-Mach number (Eq. 65) 

M — Overturning moment about the center of 
gravity 

M — logioe=0.43429 (approx.) 

N — Spin 

P — Performance parameter (Eq. 27 ) 

P — Siacci range : the distance measured along the 
line of departure to a point directly above 
the projectile (Eq. 54) 

Q — Drop of the projectile (Eq. 55) 

R — Reynolds number (Eq. 17) 

R — Range 

S — Space function (Eq. 60) 

T — Total time of flight 
T — Time function (Eq. 61) 

U — Upper limit of integration of Siacci functions 
V — Velocity 
W — Wind speed 
X — Total range of a trajectory 
Y — Altitude of an airplane above sea level 



Z — Linear deflection 
Z — Zenith angle 
i — Velocity of sound 

Diameter of the largest cross section of the 
projectile (in Gavre drag law, Eq. 23) 
c' — Damping coefficient due to yawing moment 
and crosswind force (Eq. 51) 
c" — Damping coefficient due to drag (Eq. 33) 
d — Caliber, or maximum diameter of projectile 

f — Drag coefficient (in the Gavre drag law, F.q. 

'23) 

g — Gravitational acceleration 
h — Logarithmic rate of decrease of air density 
with altitude : 0.000,031,58 per ft 
h — Yawing moment damping factor 
i — Form factor (Eq. 21) 

k — Logarithmic rate of decrease of sonic velocity 
with altitude 0.000,003,01 per ft 

k — Ratio of American to French drag function 
(Eq. 24) 

ni — Mass of the projectile 
n» — Slope of the tangent to the trajectory 
n — lx>garithmic rate of change of muzzle velocity 
with projectile weight 
p — Space function (Eq. 28) 
p — Projectile weight 
qi — Inclination function (Eq. 30) 
q — Altitude function (Eq. 31) 
s — Stability factor 

t — Time of flight (sometimes the symbol t* is used 
to distinguish the time of flight from the 
Siacci time function) 
t — Time function (Eq. 29) 
u — Velocity of the projectile relative to the air 
u — Component along the line of departure of the 
velocity relative to the air : the argument of 
the Siacci functions 

v — Velocity of the projectile relative to an inertial 
system 




w — True air speed of the airplane in aircraft fire 
x — Horizontal distance from the origin of the 
trajectory 

y — Vertical distance from the origin of the tra- 
jectory 
z — Drift 

A (capital delta) or A' — Increment of (followed 
by a variable) 

A — Density of air (in Gavre drag law, Eq. 23) 

7 (gamma) — Factor inversely p roportional to the 
ballistic coefficient (Eq. 66 ) 

9 (delta) — Yaw 

« (epsilon) — Angle of site (Eq. 72) 

A (theta) — Angle of inclination of the trajectory 
x (kappa) — Crosswind force damping factor 
|i (mu) — Viscosity of the air 
p (rho) — Density of the air 
x (tau) — Temperature of the air 
f (phi) — Angle of departure 
w ( omega ) —Angle of fall 
Superscripts 

One dot (example: x) — Derivative with respect 
to time. 

Two dots (example: y) — Second derivative with 
respect to time. 

Subscripts 

0 — At t = 0, y = 0, or 8 = 0. 



1 — Pertaining to, or relative to. Projectile Type 1 

2 — Pertaining to, or relative to, Projectile Type 2 
S— Pertaining to, or relative to, Projectile Type 8 
.1 — First revision 

.2 — Second revision 

C — Due to 1% increase in ballistic coefficient 
D — Drag 

H — Due to height of target 
R — Range component of 
e — Effective 

p — Due to projectile weight 
s — Summital : at the summit of a trajectory, where 

y = 0 

s — Standard : pertaining to standard conditions 
t — Pertaining to, or relative to, some Typical Pro- 
jectile 

v — Due to 1 unit increase in muzzle velocity 
w — Due to wind 
z — Cross component of 
9 — Yaw 

t — Due to increase in the temperature of the 
atmosphere 

? — Due to an increase in angle of departure 
— Terminal : where y = 0 on the descending 
branch of the trajectory 



v 



TRAJECTORIES, DIFFERENTIAL EFFECTS, AND DATA 

FOR PROJECTILES 



1. Introduction to Extxuok Ballistics 

Ballistics is the science of the motion and 
effects of projectiles. Exterior Ballistics deals 
with that part of the motion between the pro- 
jector and the point of burst or impact. 

The projector may be any sort of weapon: a 
sling-shot, a bow, a mortar, a rifle, a rocket 
launcher, an airplane, etc. A projectile may be 
any body that is projected by such a weapon. 
However, in this pamphlet, the term "projectile” 
will be restricted' to those missiles that are in free 
flight through the air after being propelled by a 
gun or a rocket motor. This term does not in- 
clude bombs, which are released from airplanes; 
they will be considered in another pamphlet. 

The motion of a missile in flight is really 
quite complicated. Even a fin-stabilized pro- 
jectile or bomb that is not spinning about its 
longitudinal axis is swinging about a transverse 
axis, and this affects its motion. The nose of a 
spinning shell precesses about the tangent to its 
trajectory, and its attitude gives rise to forces that 
not only affect its motion in a vertical plane but 
also move it out of this plane. To simplify the 
calculations, however, it is usually assumed that 
a missile acts like a particle subject to two forces : 
(1) gravity, which is directed vertically down- 
ward, and (2) drag, which is directed opposite to 
the direction of motion. Gravity is assumed con- 
stant : until 1956, the gravitational acceleration was 
usually taken to be 32.152 ft/sec* (9.80 m/sec*) ; 
and was used for all tables re f err e d to in this text ; 
now, the standard value is 32.174 ft /sec* (9.80665 
m/sec*). 1 The drag depends on the size of the 
missile, its velocity, the density of air, and other 
factors. Under these conditions, the motion fol- 
lows a smooth, planar path, which is called an 
"exterior ballistic trajectory,” or “trajectory” for 
short. 

The initial conditions for the trajectory of a 
particular missile are its initial velocity and angle 
of departure.* For a projectile, die ini tial velocity 

•The term angle of departure and the term elevation 
as used in this pamphlet repres e nt the ancle from (he 
koeiaontal to the line of departure and therefore would be 
m ore acc urately and completely terme d quadrant angle of 
departure . 



is the muzzle velocity and angle of depart i . < is 
approximately the elevation of the weapon 
(fictitious values of these quantities may be used 
in the case of a rocket that burns outside the 
launcher). A bomb is released with the speed of 
the airplane, generally in a horizontal direction. 
If a gun is mounted in an airplane or other 
moving vehicle, the initial conditions arc r' :cted 
by the motion of the vehicle. In any case, it is 
customary to compute trajectories for a few se- 
lected values of initial velocity and elevation. 

Certain standard conditions are specified for a 
normal trajectory : not only the initial velocity and 
angle of departure, but also the weight of the 
projectile and a factor that depends on its shape ; 
the density, temperature, and lack of motion of 
the air, and the height of the target relative to the 
projector. In order to determine the effects of 
variations from these conditions, it may be con- 
venient in some cases to compute non-standard 
trajectories; however, they can usually be esti- 
mated in other ways, as will be explained later. 

2. Vacuum Trajectories 
a. Equations of Motion 

A rough idea of an exterior ballistic trajectory 
may be gained by neglecting the air resistance, 
considering only the effects of gravity. Then if 

x is the horizontal distance from the origin, 
y is the vertical distance from the origin, 
t is the time of flight, 
g is the gravitational acceleration, 

and dots denote derivatives with respect to time, 
the differential equations of motion are 

x=0, y = — g. (1) 

The initial conditions, denoted by the subscript 0, 
arc 

x. = v.cos $*,y» = v„ sin •&>, x, = 0, y D = 0, 

where v is the velocity, 

S the angle of inclination. 

By integrating equations (l), we find that, at 
any time, 

x = v. cos y = v„ sin — gt, (2; 

and by integrating (2) , 
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x =r (v 0 cos -So)t, y = (vo sin $o)t — gt 2 /2. (3) 
b. Characteristics of Trajectories 

From Equations (3), we find that the time at 
any horizontaJ range is 

t = x/v, cos $<> (4) 

and hence the ordinate is 
y = x tan -So — gx 2 / 2v 0 2 cos 2 $„• (5) 

This is a parabola with a vertical axis. 

The total range X is found by letting y =0 in 
(5) and solving for x: 

X = (2v c 2 /g) sin <S 0 cos = (v„ 2 /g) sin 2^ 0 . (6) 
This shows that the range is a maximum when 
•So = 45°. 

The total time of flight T is found by substi- 
tuting (6) for x in (4) : 

T = (2v*/g) sin -So. (7) 

The components of the terminal velocity are 
found by substituting (7) for t in (2): 

X = v 0 cos $o» Y = — Vo sin * (8) 

Hence the terminal velocity is 
v„=Vo, (9) 

and the angle of fall (considered positive) is 
« = • So. (10) 

By differentiating equation (5) with respect to 
x, we find that the slope of the trajectory at any 
point is 

^ = tan <S 0 — gx/v 0 * cos* •$*. (11) 

At the summit, the slope is 0 ; hence the horizontal 
range to the summit is 

x, = (v,*/g) sin <S 0 cos $« = X/2. (12) 

By substituting this for x in (5), we find the max- 
imum ordinate 

y, = (v 0 */2g) sin 2 -So = ( 1 — COS 2$ c ) v,*/4g 

= (g/8) T*. (13) 

By substituting (12) for x in (4), we find the 
time to the summit 

t* = (vo/g) sin = T/2. (14) 

3. Air Resistance 

a. Forces Acting on a Missile Moving m Air 
A missile moving in air encounters other forces 
besides gravity. The principal components of the 
force produced by motion through air are drag, 
crosswind force, and, for a spinning projectile, 
Magnus force. 

( 1 ) Drag is the component in the direction op- 



posite to that of the motion of the center of 
gravity. This is the only component acting on a 
non-spinning projectile that is trailing perfectly. 

(2) Crosswind Force in the component per- 
pendicular to the direction of motion of the center 
of gravity in the plane of yaw. The yaw is the 
angle between the direction of motion of the 
center of gravity and the longitudinal axis of the 
missile. This component is the cause of the drift 
of a spinning shell. 

(3) Magnus Force acts in a direction per- 
pendicular to the plane of yaw. This is the force 
that makes a spinning golf ball or baseball swerve. 

b. Simplified Assumptions 

The motion caused by these forces is explained 
in detail in the Ordnance Corps Pamphlet, “De- 
sign for Control of Flight Characteristics.”* In 
general, the motion of the center of gravity is 
affected by the motion about the center of gravity 
and occurs in three dimensions. To simplify the 
calculations, it is usually assumed that the only 
forces acting on the missile in free flight are 
drag and gravity. The crosswind force is implied 
in the drift, but this is determined experimentally 
as a variation from the plane trajectory. This 
approximation is satisfactory except at extremely 
high elevations, where the summital yaw is very 
large. 

c. Drag Coefficient 

The drag encountered by a projectile moving in 
air is defined by the equation of motion 
dv 

m = — D — mg sin •$, (IS) 

where 

m — is the mass of the projectile, 
v — the velocity of the projectile relative to an 
inertial system, 
t — the time, 

D — the drag, 

g — the gravitational acceleration, 

$ — the angle of inclination. 

The drag coefficient is defined by the relation 
Kd = D/pd*u>, (16) 

where 

Kx> — is the drag coefficient, 
p — the air density, 
d — the caliber or maximum diameter, 
u — the velocity of the projectile relative to 
the air. 

The drag coefficient is a function of the Reynolds 
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number R and the Mach numlier M. The 
Reynold# numlier is defined as 

R = udp/p (17) 

where |i is the viscosity of the air, and the Mach 
number as 

M = u/a, (18) 

where a is the velocity of sound in the surround- 
ing medium.* The viscosity is related to the 
friction of the air ; the velocity of sound, to its 
elasticity. At extremely low velocities, the 
Reynolds number causes greater variation, but 
the effect of the Mach number is larger when it 
exceeds 0.5. The angle between the direction of 
motion of the center of gravity and the axis of 
symmetry of the shell, called yaw, also affects the 
drag coefficient, but this is usually stripped out of 
the resistance firing data. 

Figure 2 (p. 22) is a graph of drag coefficient 
vs Mach number for the four typical projectiles 
whose shapes are shown in Figure 1 (p. 21). 
These projectiles are described and the basic 
firings explained in Section 4. It should be noted 
that subscripts are used to distinguish the typical 
projectiles to which the drag coefficients pertain ; 
a second subscript is sometimes added after a 
decimal point to denote a revision. It is evident 
that all the drag coefficients increase rapidly in 
the vicinity of Mach number 1. The main differ- 
ence between curves is in the ratio of the maxi- 
mum, which occurs between Mach numbers 1.0 
and 1.5, to the minimum subsonic value. Although 
the drag coefficient has been determined and tab- 
ulated for several other projectiles, these four 
represent the principal shapes for artillery pro- 
jectiles. 

Kd is a non-dimensional coefficient, and a con- 
sistent set of units should be used in Equation 
(16). However, it is customary to express the 
density as a ratio, the caliber in inches, and the 
velocity in feet per second ; then, since the drag 
is in poundals, the drag coefficient is expressed in 
pounds per square inch per foot and is denoted by 
B. Until 1956, the standard air density was 
taken as 

p„ = 0.07513 lb/ft 3 = 5.217 X 10~ 4 lb/in 2 -ft * 

The practical drag coefficient, 

B = 5.217 X 10- 4 K d , (19) 

for Projectile Types 1, 2 and 8 is tabulated as a 
function of velocity for the standard velocity of 
sound, 

ap= 1120.27 fps.* 

* In 1956, the Department of Defense adopted the stand- 
ard atmosphere of the International Civil Aviation Organ- 
ization. 1 In this atmosphere, at the surface, the density is 
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The non-dimensional drag coefficient Kd for the 
90-mm HEAT Shell T 108 is tabulated as a func- 
tion of Mach number. (See References 3, 4, 5 and 
6, which may be obtained from the Ballistic Re- 
search Laboratories, Aberdeen Proving Ground, 
Maryland). 

d. Drag Function 

The drag function G is computed by the formula 
G = Bu (20) 

with u in feet per second. For convenience in 
computing trajectories, this has been tabulated 
as a function of u 2 /100 with u in meters per 
second and either feet or yards per second, for 
Projectile Types 1, 2 and 8. (See References 7, 
8 and 9.) 

e. Form Factor 

It is not practical to determine the drag co- 
efficient as a function of Mach number for all 
projectiles. Instead, the drag coefficient of a pro- 
jectile whose shape differs only slightly from one 
of the typical projectiles is assumed to be pro- 
portional to the typical drag coefficient. If Kn t 
is the drag coefficient of a Typical Projectile, the 
form factor of a projectile whose drag coefficient 
is K d is the ratio 

i t =K D /K Dt . (21) 

The form factor may be determined from the 
time of flight at short ranges, or from the range 
to impact or the position of burst at a few eleva- 
tions. If it is desired to estimate the form factor of 
a projectile without firing it, the Ordnance Corps 
Pamphlet, “Design for Control of Flight Char- 
acteristics” 2 should be consulted. It explains the 
effects on drag of variations in length and radius 
of head, meplat diameter, base area, length of 
shell, and yaw. 

Generally, the form factor of a projectile with 
a square base and an ogive less than 1.75 calibers 
high should be referred to Projectile Type 1 ; that 
of a projectile with a square base and an ogive 
more than 1.75 calibers high, to Projectile Type 
8; that of a projectile with a boattail, to Projectile 
Type 2; and that of a finned projectile, to the 
T108 Shell. However, there are exceptions to 
these rules :e.g., K D , is used for fin-stabilized mor- 
tar shells at low velocities, and Koa fits the calcu- 
lated drag coefficient for some longer fin-stabilized 
shells. The drag coefficients of bombs are treated 
elsewhere. 10 

0.076,475 lb/ft* and the velocity of sound is 1116.89 fps; 
also the viscosity is 1.205 x 10~' lb/ft-sec, but the effect 
of Reynolds number is not considered in computing tra- 
jectories. 
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}. Ballistic Coefficient 

The ballistic coefficient relative to the tabulated 
drag coefficient is defined by the formula 

C t = m/i,d 2 , (22) 

where m is the mass of the projectile in pounds, 
d is the caliber (or maximum diameter) 
in inches. 
g. Air Density 

The air density is a function of altitude. The 
exponential function adopted by the Ordnance 
Corps* is tabulated in References 7 and 11. The 
ratio of the density at an altitude y to that at the 
surface is H(y). 

4. Description of Typical Projectiles 

This pamphlet contains data that pertain spe- 
cifically to four types of projectiles: Types 1, 2 
and 8 and the 90-nun HEAT Shell T108. They 
can be applied to other shapes by using a form 
factor, as explained in paragraph 3e. 

a. Projectile Type 1 

Most of the service projectiles used by the 
European artillery during the last third of the 
nineteenth century had square bases and ogival 
heads of various apical angles. In Russia, during 
the years 1868 and 1869, General Magenski con- 
ducted some experiments at St. Petersburg to de- 
termine the law of air resistance of projectiles with 
a semi-apical angle of 48°22' (1.49-caliber radius 
head) at velocities from 172 to 409 m/s. In Eng- 
land, from 1866 to 1880, Bashforth conducted sim- 
ilar experiments with projectiles of the same form 
fit may be that the tables were reduced to this 
form) at velocities from 148 to 825 m/s. In Hol- 
land, Colonel Hojel conducted a few such experi- 
ments during 1883; the exact form is uncertain, 
but the ogive was evidently somewhat longer than 
those of the Russian experiments. In Germany, 
the Krupp factory conducted a large number of 
air resistance experiments at Meppen Proving 
Ground, starting in 1881 (the principal results 
were reported by Siacci in the “Rivista” for 
March 1896) ; the velocities varied from 368 to 
910 m/s, and three different types of projectiles 

* The tabulated trajectory data based on the old stand- 
ard atmosphere give results that are nearly the same as 
those based on the new atmosphere providing a suitable 
change is made in the ballistic coefficient. For example, a 
280-mm shell, fired with a muzzle velocity of 1795 fps, at 
an elevation of 45°, attains a range of 15,790 yds.; the 
ballistic coefficient is 3.62 with ICAO atmosphere and 
3.50 with the old Ordnance Corps atmosphere. This func- 
tion was used in computing trajectories before 1956, when 
the Department of Defense adopted the standard atmos- 
phere of the International Civil Aviation Organization. 



were used: one was the same as the Russian (or 
the results reduced to the same), but the others 
probably had longer ogives. In France, a com- 
mission of the naval artillery conducted many such 
experiments at Gavre Proving Ground from 1873 
to 1898 at velocities from 118 to 1163 m/s; they 
fired ogival projectiles with three semi-apical 
angles: 45°55' (1.64 CRH), 41°40' (1.98 CRH), 
and 31°45' (3.34 CRH). For further details, see 
References 12 to 16. 

In analyzing the results of the experiments, 17 
the Gavre Commission of Experiments assumed 
that the drag could be expressed in the form 

D = — a 2 V l f(V), (23) 

where 

D — is the drag (kg), 

A — the density of the air (kg/m 8 ), 

g — the gravitational acceleration (m/sec 2 ), 

a — the diameter of the largest cross section of 
the projectile (m), 

V — the velocity (m/s). 

The drag was treated as a function of velocity 
without taking account of the velocity of sound, 
although the dependence of drag on the velocity of 
sound was strongly suspected: the coefficient f 
was plotted as a function of V alone. After 
studying the results obtained with projectiles of 
various ogival angles, the Commission confirmed 
the law deduced by Helie in 1888 that f(V) is 
proportional to the sine of the ogival angle, at 
least for values from 40° to 90°, within the experi- 
mental errors. 

In 1917, Chief Engineer Gamier of the Gavre 
Commission published a table 18 of log B and the 
corresponding logarithmic derivative, with ve- 
locity as argument: this B is equivalent to f(V) 
for an ogival projectile with a 2-caliber radius 
head. This ideal shape is what we call Projectile 
Type 1 ( see Figure 1 ) . The Ordnance Depart- 
ment of the United States Army then prepared a 
table of the logarithm of the drag function 19 

G = kBV, (24) 

with V expressed in meters per second, and the 
factor k to take account of differences in defini- 
tions and standards. Below 600 m/s, k = 0.001 14 ; 
but at higher velocities, it is variable. G was 
called the Gavre function; later, the symbol was 
changed to Gy to denote that it pertains to Pro- 
jectile Type 1. 

In the table of log G, the argument, V 2 /100, 
varies from to 0 to 33,000 (equivalent to 1817 
m/s or 5961 fps). G has been tabulated for 
V 2 /100 up to 25,000 (1581 m/s or 5187 fps). 
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The product Gv has also been tabulated as a 
function of v in feet per second up to 6000 fps. 10 

b. Projectile Type 2 

A sketch of Projectile Type 2 is shown in 
Figure i. This projectile has a 6° boattail y% 
caliber long and an ogivo -conical head 2.7 cali- 
bers long. The experimental 4.7-inch Shell El 
had this shape. However, the resistance firings 
were actually conducted with two types of 3.3-inch 
shell: Type 155 with a 5° boattail and Type 157 
with a 7 C boattail. The results were so close that 
they were grouped together to represent one type. 

These projectiles were fired from 1922 to 1925 
at velocities from 655 to 3187 fps. The results of 
these firings were fitted by the tabulated drag 
function G* (this is defined in paragraph 3d). 
Later, it was extrapolated to 10,000 fps by means 
of theoretical considerations. For the sake of 
greater smoothness and accuracy, two revisions 
have been made: G*.i, tabulated to 6000 fps, and 
Gj, 2 . to 7000 fps. Gs.i is approximately the same 
as G*, but Gs .2 is lower at velocities above 2000 
fps. 

c. Projectile Type 8 

A sketch of Projectile Type 8 is shown in Fig- 
ure 1. Tnis projectile has a square base and a 
secant ogive 2.18 calibers high. Its ogival radius 
is 10 calibers, which is twice that of a tangent 
ogive of the same hdg.it. 

The drag function G, was calculated from data 
obtained oy the British Ordnance Board from 
inree sources: 

( 1 ) Aieasuremems of a model in the National 
Pnysica. Laboratory wind tunnel at Mach numoers 
from 0.22 to 0 .if*. 

(2 Resistance firings of 2-pounder sod: at 
Macii numoers from 0.65 to 2.8. 

•'3) Resistance firings of o-inch shell at Macr. 
numbers from 0.55 to 4 l. ; . 

G» was first tabulated to 5000 fps. Later, a 
revised table, caned G*.i, was extended to 6000 
fps. Toe corresponding drag coefficients are ap- 
proximately proportional to the British “1940 
’-aw,” which was first tabulated to 5000 fps and 
later revised above 3750 fps and extended to 6000 
fps by a power law. 

d. HEAT Shell, 90-mm, T108 

A sketch of the 90-mm High Explosive Anti- 
tank Shell T108 is shown in Figure 1. This pro- 
jectile has fins one caliber in span at the end of a 
long boom and a conical ogive about 2 calibers 

high. 

The drag of this projectile was determined at 
velocities from 700 to 2700 fps by free flight firings 



in the large Spark Range of the Exterior Bal- 
listic Laboratory at Aberdeen Proving Ground 
The drag coefficient (defined in paragraph 3d 
was fitted with a series of analytical expressions, 
which were used in tabulating it as a function o- 
Mach number from 0.1 to 2.7. 



5. Trajectories in Air 
a. Equations oj Moucn 

Under standard conditions, the equations c 
tion are 



GH . ' 

x = — ~-x, y = 



GH . 



where x is the horizontal range, 
y — the altitude, 

H — the air density ratio, 
g — the gravitational acceleration (32.i:2 
ft /sec 2 ), 

and dots denote derivatives with respect to time. 
The initial conditions are 

Xo = 0, y 0 = 0, Xo = v 0 cot b,„ y P = v 0 sin 

where v 0 is the initial velocity, 

£<, — the angle of departure. 

These equations cannot be solved analytical'^, 
but they can be solved by a method of numerical 
integration explained in References 21, 22, 23 and 
24. 



b. Characteristics of Trajectories 

(1) Ballistic Tobies. Trajectory data have 
been tabulated tor Projectile Types 1 &t;d 2 and 
for some bomb shapes. 

The “Exterior Ballistic Tables Based on Nu- 
merical integration’* (Reierences 7, 25. and 26' 
pertain to Projectile Type I. The merer is tne 
unit of lengtiv Tne trajectories were computed 
forward and backward from the summi.: and cut 
off at several altitude.-. It C* is the ‘‘rummita: 
ballistic coefficient” of . trajectory and y, the max- 
imum ordinate of a partial trajectory, the ballistic 
coefficient C of this partial trajectory satisfies .We 
relation 



logioC = logicC, — 0.000,045 y,. (26) 

For each y„ C,, and summital velocity v„ 
Volume I gives the horizontal range, t;me. and 
velocity components at the beginning and end 
of the partial trajectory. If the muzzle velocity, 
angle of departure, and ballistic coefficient are 
known, Volume II may be used to find the 
time of flight to the summit, t,, the velocity at the 
summit, the horizontal range to the summit, x*, 
and the maximum ordinate. Then C. can be com- 
puted and trajectory data can be found in Volume 
I for the given C, and v f at various altitudes. 




However, terminal values can be found more con- 
veniently from Volume III, which gives the hori- 
zontal range X. time of flight T, angle of fall 
u>, and terminal velocity v w ; the arguments are 
the same as in Volume II: the muzzle velocity 
vanes from 80 to 920 meters per second, the 
angle of departure from 0° to 80° (at intervals of 
1° from 0° to 30°. and 2 ° from 30° to 80°). the 
liallistic coefficient from 1.000 to 17.78 (at inter- 
\ a Is of 0.05 in log C ) . 

The “Brief Exterior Ballistic Tables for Pro- 
lertile Yvpe i” (Reference 27) contains values 
that were obtained by interpolation in Volume 
ill and conversion to English units (yards and 
feet per serond i . The muzzle velocity varies from 
2fX) to 3000 feet per second in 100 foot per second 
increments. The angle of departure varies from 
15° to 75° in 15° increments. Obviously, inter- 
polation for angle of departvire in this table is 
difficult ; nevertheless, the table is useful for many 



purposes. The table of range for an elevation of 
45°. which is nearly’ the maximum range, is repro- 
duced here ( Table 1 ) . 

The "Exter.or Ballistic Table for Projectile 
Type 2” (Reference 28) consists of three parts: 
Table I, Summital and Terminal Values; Table 
II. Coordinates as Functions of Time of Flight; 
and a large graph of Maximum Range vs. Ballistic 
Coefficient. The tables have been extended to 
include data for every 15° from 0° to 90°. T r.e 
muzzle velocity varies from 12C0 to 5200 feet ; . 
second, and the ballistic coefficient from 0.50 to 
10.00. The coordinates are given in yards, the 
time in seconds, and the velocity in feet per second. 
Table I gives the time to the summit, the hori- 
zontal range to the summit, the maximum 
ordinate, the time of flight, the total range, the 
angle of fall, and the terminal velocity. Table il 
gives the horizontal range and altitude at every 
5 seconds of time. Unfortunately, the velocity 



TABLE I 

MAXIMUM RANGE FOR PROJECTILE TYPE 1 
ELEVATION 45° 

Range in yards. Muzzle Velocity (v„) in fps 



I-Og Oi 

c, 

v o 


0.000 


0.050 


0.100 


0.150 


m 






EE9I 


m 


RH 




0.550 


0.600 1 


m\ 


1.000 


1.122 


1 .259 


1.413 


1 .585 


1.778 


1.995 


2.239 


2.512 


2.818 


3.162 


3.M8 


3.981 

--- 


4.467 


200 ; 


307 


398 


399 


400 




403 


404 


405 


406 


407 


408 


408 


jWfjl 


409 


300 


K47 


856 


864 


871 


877 


883 


888 


892 


897 


000 


004 


907 




913 


400 ! 


1410 


1441 


1462 


1481 


1499 


1514 


1528 


1541 


1553 


1564 


1573 


1582 




Esa 


500 


2077 


2121 


2162 


2201 


2237 


2270 




2328 


2354 


2378 


2399 


2418 


2436 


Cal 


too 


27 M6 


2862 


2934 






3126 


3181 


3232 


3279 


3322 


3361 


3397 


3428 


3461 


700 


3521 


:Xi38 


3740 


3856 


3958 


4053 


KlLifl 


4222 




4372 


4438 


4500 


4556 


4606 


MOO 


4247 


4411 


-'■•570 


4724 




Kj£| 


5141 


526.7 


5384 


5494 


5595 


5690 


5777 


5859 


f >on 


4020 


5130 


5340 


-5556 


5755 


5947 


6133 


6310 


647C 


6634 


6783 


6923 


7054 


7174 




54MH 


5757 


6024 


6286 


6542 


6792 


7037 


7272 


7498 


7714 


7919 


8115 


8299 


8472 


MTS'S 


5026 


6238 


1.553 


6866 


7176 


7482 


7784 


8070 


8366 


8644 


8913 


9172 




9658 


1200 il 6251 


6603 


60.58 


7314 


7671 


8027 


8381 


8733 




9417 


9749 






Eml 




6521 


6005 


7293 


7689 


8087 


8488 


8801 


9292 


9603 






10872 


11253 


11626 




6759 


7172 


7504 


8025 


8460 


8905 


9353 


9804 






11165 


11616 




WWT : 11 


1500 1 


6979 


7419 


7872 


8336 


8810 


9294 


9787 








11814 




12845 


13360 


uoo 


7185 


7051 


8134 


8630 


9141 


0664 


10200 


10747 


11302 


11860 


12442 


13021 


13605 


14194 


1700 j 


7>1 


7873 


8384 


8012 


9458 


10021 


10500 


11193 


11798 


12421 


13056 


■kriijl 


14356 


Kml 


IMHO ; 


7569 


8088 


8627 


9187 


9767 


10368 


10080 


11630 


12288 


12967 


13662 


14374 


15099 


15839 


1000 ' 


7752 


8206 


8863 


04.54 


10068 


10707 


11371 


12059 


12771 


13505 


14262 




15839 


16659 | 


2UI) ! 


7031 


8408 


0003 


9715 


10364 


11042 


11740 


12483 


13247 


14037 


14850 








2100 ',i 81(16 


8608 


9320 


9972 


10655 


11372 


12121 


12903 


13719 


14567 


15454 


16372 


17323 




2200 i; K27X 


8804 


9542 


10225 


10942 


11697 


12489 


13310 


14188 


15096 


16040 


17038 




19136 


2300 ! 


8447 


0086 


0760 


10473 


11225 


12018 


12853 


13732 


14655 


15624 


16644 


17706 




19975 


I) ! 


8612 


0274 


0075 


10718 


11.505 


12336 


13215 


14143 


15121 


16152 


17239 


18375 


19574 




2500 


8774 


9450 


10188 


10961 


11782 


12652 


13575 


14552 


15585 


16679 


17834 


19048 




I 21675 | 




8935 


9642 




11201 


12056 


12065 


13931 


14958 


16048 


17205 


18430 


19725 




22539 




Era 


9823 




11438 


13237 


13275 


14285 


15362 


16510 


17731 


19027 


20405 






. 2WO 


! 9247 






11672 


12595 


13582 


14637 


15764 


16970 


18257 


19625 


21091 






2! if 10 1 


9309 






11903 


12861 


13886 


14087 


16165 


17429 


18782 


20225 


21778 


23428 


25184 


3000 


9549 

j 


10351 




12132 


13123 


14187 1 15334 


16563 


r.s85 


10303 


20823 


22464 


24213 
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TABLE I (Continued) 

MAXIMUM RANGE FOR PROJECTILE TYPE 1 
ELEVATION 45° 



Range in yards. Muzzle Velocity (v,) in fps 




MAXIMUM RANGE OF PROJECTILE TYPE 2 
Range in Units of 10 Yards 




components are not tabulated. The graph has a 
curve of maximum range vs. ballistic coefficient 
for every 400 feet per second muzzle velocity ; it 
also lists the angle of departure for maximum 
range. A table f* the maximum range is included 



in this pamphlet (Table II). 

Some trajectories have been computed for Pro- 
jectile Type 8 at an angle of departure of 45°. 
The range, which is nearly the maximum range, 
is tabulated here (Table III). 
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TABLE 111 

MAXIMUM RANGE FOR PROJECTILE TYPE 8 
ELEVATION 45° 



Range in Units of 10 Yards 



Muule Velocity tpe 


Log C» 


-0.2 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


j 1200 


521 


665 


815 


061 


1091 


1201 


1288 ' 


1 2000 


680 


023 


1223 


1581 


1034 


2403 


2826 1 


2800 


822 


1168 


1650 


2287 


3124 


4231 


5251 


3900 


064 


1432 


2145 


3162 


4865 


6936 


8741 


4400 


1108 


1720 


2726 


4470 


7492 


10884 


13661 t 


5200 


1258 


2040 


3424 


6604 


11603 


16582 


20310 

l 



When bombs were dropped at low speeds, tra- 
jectory data based on the Gavre drag function Gi 
were accurate enough for calibrating the sights. 
In fact, the descending branches of the trajectories 
in Volume I of the “Exterior Ballistic Tables" 
could be used for this purpose. Actually, some 
“Bomb Ballistic Reduction Tables” were pub- 
lished, which were more convenient for this pur- 
pose. However, at transonic speeds, the drag 
functions of modem bombs differ considerably 
from Gj. Drag coefficients have been determined 
for several shapes, but three of them have been 
selected as bases of "Bomb Ballistic Tables" (Ref- 
erence 10) : drag coefficient A was determined 
from wind tunnel tests of a shape designed by the 
National Advisory Committee for Aeronautics, and 
is probably the lowest drag coefficient that any 
bomb will have ; drag coefficient B was determined 
from range bombing of the 3000-pound Demoli- 
tion Bomb Ml 19, and is suitable for other mem- 
bers of the series including the 750-pound Bomb 
M118 and the 10,000-pound Bomb M121; drag 
coefficient C was determined from resistance fir- 
ings of the 75-mm Slug Mark I, and is probably 
the highest coefficient that any bomb will have. 

The arguments of the “Bomb Ballistic Tables” 
are release altitude ( 1000 to 80,000 feet above sea 
level), true air speed of release in level flight (250 
to 2000 knots), and performance parameter (0 
to 1.5). The performance parameter P is inverse- 
ly proportional to the ballistic coefficient; it is de- 
fined by the formula 

P = lOOd’/m, (27) 

where d is the maximum diameter of the bomb in 
feet. 



m — the mass of the bomb in pounds. 

The tables give the horizontal range from release 
to impact (feet), the time from release to impact 
(seconds), the striking velocity (ft/sec), and the 
striking angle (degrees). The drag coefficients 
are tabulated as functions of the Mach number in 
the appendix ; they are graphed in the intro- 
duction. 

(2) Graphs. Figures 3, 4 and 5 are graphs of 
the total time of flight for an elevation of 45°. The 
time of flight is useful in estimating wind effects. 

For antiaircraft firing and for firing to the 
ground at different altitudes, points along the tra- 
jectory are needed. Figures 6 and 7 show the 
trajectories computed with G 2 and Gg at angles of 
departure of 45° and 60“ (there is a different 
graph for each muzzle velocity, each elevation, 
and each Projectile Type). Trajectories based 
on Gi are not included because no antiaircraft pro- 
jectiles have a blunt head. The T108 HEAT 
Shell is usually used only in direct fire at short 
ranges. Trajectory data for other elevations can 
be estimated from those at 45° or 60° by assuming 
that the slant range is the same for a given time of 
flight. 

Unlike the vacuum trajectories, in air, the 
horizontal range to the summit is more than half 
the total range, and the time to the summit is less 
than half the total time; however, the maximtup 
ordinate is approximately (g/8)T*. In air, the 
terminal velocity is less than the initial velocity 
and the angle of fall is greater than the angle of 
departure. At a given elevation, the total r$nge 
increases less rapidly than the square of the 
muzzle velocity; at a given muzzle velocity, the 



8 




range is a maximum at an angle somewhat differ* 
ent from 45* : lower for small ballistic coefficients, 
higher for large ones. However, the range at 45° 
is nearly as much as the maximum. 

6. Siacci Tables 

a. Spin- stabilised Projectiles 

(1) Tabular Functions. In 1880, the Italian 
Col. Siacci introduced the Space, Time, Altitude, 
and Inclination Functions, (Reference 29), which 
simplify the calculation of trajectory data under 
certain restrictions. In particular, their use will 
be explained for three conditions: ground fire at 
low superelevations, antiaircraft fire (where the 
superelevation is small), and aircraft fire. These 
and some auxiliary functions based on Gi, Gs.i 
(the second revision of Gs) and G«.i (the first 
revision of G#) have been tabulated for the Air 
Force (References 30, 31 and 32). Since these 
Siacci tables were arranged primarily for use in 
computing aircraft firing tables, they were called 
Aircraft Tables; however, they can be used for 
other purposes, as explained below, in lieu of other 
forms. In these Aircraft Tables the argument is 
the component u along the line of departure of the 
velocity relative to the air, which is assumed mo- 
tionless. The following functions and some of 
their derivatives with respect to u are tabulated : 



Altitude* 



r u q,du _ 

q Ju 



(31) 



Space 




fU 


du 


(28) a/i 


p-j 


L 


G(u) * 


Time 


t= J 


fU 


du 


(29) 

If 


l 


uG(u) 


Inclination* 


q,= j 


fV 


g du 


(30) > 


l 


u a G(u) ’ 






X — 


(Cpo/p) cos $ 


. [p(ua«/a) 



Velocity increment A'u = -^pG(u), (32) 

(K d $= 16.4 per rad*), 

Damping coefficient c" = G(u) /2u. ( 33) 

The upper limit of integration, U, is oC '0 fps. tor 
Projectile Type 1 and 7000 fps for Types 2 and 
8. K D g is the yaw-drag coefficient, defined by 
the formula 

K»=Kn B (l + K D a a ), (X) 

where Ko„ is the drag coefficient for 0 yaw, 
i — is the yaw in radians. 

(2) Ground Fire at Low Elevations 

(a) Formulas. Although the Aircraft Tables 
are not accurate in computing complete trajectories 
at high elevations, they can be used to obtain tra- 
jectory data at superelevations below 15° (for the 
derivation and application of the following for- 
mulas, see Reference 33). Here, we shall take 

p/po — as the ratio of air density to normal, as- 
sumed constant for each problem, 

the ration of the velocity of sound to 
normal, assumed constant for each 
problem. 



a/ao = \/(459 -f F) /518. 
Then the horizontal range is 



the time of flight is 

t f = (Cp«a«/pa) [t(ua„/a) -t(uA/a)J, 

the altitude is 

y = X tan — (Cpoa,/pa)» (q(ua*/a) — q(u«a*/a) ] + x(Cp,a,/pa) sec [q x (u«a«/a) ], 
and the siope of the trajectory is 

tan $ = tan — (Cpo*.*/?* 1 ) sec b a [qi (ua„/a) — q, (u^/a)]. 

The angle of departure for y = 0, may be found by formula 

sin * t = Cp^r qLu a^ a) - q (<U» . £) _ ^ 

pa Lp(uao/a) — p(u g a*/a) ' y J 

If p/p. 1 and a/a. = 1, let p = p(u), p. = p (u.), etc. 

* In most places, the altitude and inclination functions are defined aa twice these integrals. 
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( 35 ) 

( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 



Then these formulas can be expressed more 



simply : 

x — C cos (p ■— po) , (36a) 

t f = C(t_0, (37a) 

' y = x tan — C*(q — q*) + 

C sec $ 0 x q JO ) (38a) 

tan & — tan — C sec $ c (q, — q jc ), (39a) 

tnd, for y ~ 0, 

sin -5.^ C r^-Si-q,. 1 (40a) 

LP — Po J 



(h) Example 

Data: Tank Gun, 76-mm, M1A2; Shot, AP, 
M339. 



Caliber 

Mass 

Muzzle Velocity 
Form Factor 
Air Density Ratio 
Sound Velocity Ratio 



3.000 in. 
14.58 lb. 

3200 fps 
1.164 on Gm 

1.000 
1.000 



Problem: To find the range, time of flight, angle 
of departure, and angle of fall where 
the velocity is 2000 fps. 



Solution : From the Aircraft Table Based on Gs.i : 



u 


P 


t 


q 


9i 


2000 


23079.3 


6.5132 


489.0 


.06523 


3200 


15425.1 


3.4997 


161.6 


.02637 


Diff. 


7654.2" 


3.0135 


327.4 


.03886 



C,.! = 14.58/1.164 (3.000)*= 1.392 (Eq.22) 

sin = 1.392 ( 327.4/7654.2 - .02637) = 
.02283 (Eq. 40a) 

$0 = 23.26 mils, cos = .99974, sec. = 
1.00027 



tan$.= .02284 

1.392(1.00027) (.03886) = .05411 

tan d = —.03127 (Eq. 39a) 

£ = — 31.84 mils 

x = 1.392 (.99974) (7654.2) = 

10652 ft (Eq. 36a) 

t,= 1.392 (3.0135) = 4.19 sec (Eq. 37a) 



Note : Since u is the component of the velocity 
along the line of fire, the tangential velocity is 
v = u cos £«/cos 

In this example, v = 2000 (. 99974 )/.9995l = 
2000.5 fps. 

(3) Antiaircraft Fire 

(a) Formulas, These tables can also be 
used for the flat part of antiaircraft trajectories 
where the superelevation is less than 15* (for the 
derivation and application of the following formu- 
las, see Reference 33). Let p/p 0 and possibly C 
be known variable functions of x, but assume a 
constant value of a/a« (either 1 or an average). 
Since p/p 0 is a function of y, denoted by H, its 
relation to x can be found by the approximate sub- 
stitution 

y ~ xtan$„. (41) 

Let m = tan $, m, = tan 6». . 

Then the following trajectory data can be found 
by Simpson’s rule or some other method of nu- 
merical integration: 



x 

p (ua*/a) = p(u«a*/a) 4- sec dx, (42) 

o 



t, = sec^ 0 J* , 
o 



(43) 



m 



= m .-pi£^dx. (44) 



= J*m dx. 



( 45 ) 



Check : x tan $ 0 = 10652 (.02284) = 243.21, 

- Ch (q - q 0 ) = ( 1-392)* (327.4) = - 634.40 

C,.i secK X q, 0 = 1-392 (1.00027) (10652) (.02637) = 391.11 

y = 
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0.00 (Eq. 38a) 




(b) Example 



Dali: Automatic (AA) Qua, 40-nun, M2A1; Shot, AP, M81 
Muulo Velocity : 2870 fpa 
Bollktie CooHMaot : 0.574 on Q«., 

AagU at DiptrUn : 1200 mik 
Air Density Ratio : H-a* 

Bound V4ocity Ratio: 1.000 

Problem: To And the time of flight and altitude for horizontal range* from 0 to 2400 ft. 
Solution: (See remark* below and urn the Aircraft Table Baaed on 0,.i). 



X 

ft 

1 


Mhy 

Mt 

2 


H 

3 


fudx 

4 


P(u) 

6 


u 

6 


2613.1 

u 

7 


*f 

MC 

8 


rry 

9 


iff y* 
10 


m 

11 


y 

ft 

12 


0 




1.0000 


0.00 


17430.6 




0.0105 




0.8200 


0 


2.4142 


HK1 


300 




.0774 




















000 


.01087 


.0563 


680.40 




2446.4 


1.0681 


0.604 


1.1408 


501 


2.3062 


1443 


000 




.0337 




















1200 


■ 


.0126 


1146.76 


22800.1 


2060.6 


1.2681 


1.295 


1.6081 


1416 


2.3687 


2872 


1500 


B <Ir ■ M 


.8010 




















1300 


B v/’TiB 


.8718 


1681.06 


25006.6 


1716.3 


1.5225 


2.132 


2.3180 


2503 


2.3308 


4282 


2100 


.00063 


.8621 














- 






2400 


.07046 


.8328 


2103.26 


27424.2 


1413.7 


1.8484 


3.143 


3.4166 


4314 


2.2765 


5664 




























Remarks : 

Col. 1. The interval of x should be chosen small 
enough so that the integrations may be performed 
accurately by Simpson's rule: not more than 600 
feet. 

Col. 2. M = logio e = 0.43429. h = 0.000,031,- 
58 ft -1 . The estimate of y = x tan A 
slightly smaller factor would be better at the 
longer ranges. 

Col. 3. For convenience, the formula H = anti- 
colog 10 Mhy was used. The density function 
H(y) has been tabulated (Reference 11) ; it gives 
the same values. The ICAO “Standard Atmos- 
phere-Tables and Data” 1 gives slightly different 
values. 

Col. 4. Simpson's rule was used, but the trap- 
ezoid rule gives 2193.30 for the last value. 

Col. 5. The first value was found in the Aircraft 
Table for u*, = 2870 fps. The other values were 
computed by Eq. (42) with 

sec VC,., = 2.6131/0.574 = 4.5524. 

Col. 6. The values of u after the first were found 
in the Aircraft Table by inverse interpolation, 
using the first two terms of Taylor’s series. 

aec = 2.6131. This is multiplied by 
1000 to avoid writing zeros. 



Col. S. Eq. (43). Simpson’s rule gives 3.131 
for the last value, but the trapezoid rule was used 
in order to obtain intermediate values. 

Col. 9. The square of Col. 7. 

Col. 10. Note the factor 10* in these values. 
Simpson’s rule gives 4259 for the last value, but 
the trapezoid rule was used so as to obtain more 
values. 

Col. 11. Computed by Eq. (44) with 

nrv, = tan g = 32. 1 52 ft /sec 1 . 

Col. 12. Computed by Eq. (45). Simpson’s rule 
gives 5666 ft for the final value. 

(4) Aircraft Fire 

(a) Formulas. These tables are especially 
useful for aircraft fire (for the derivation and ap- 
plication of the following formuas, see References 
34 and 35). Before using them, the initial con- 
ditions must be determined. Let 

A — be the azimuth, measured from the direc- 
tion of motion of the airplane, 

Z — the zenith angle, measured from the upward 
vertical, 

v 0 — the muzzle velocity (relative to the gun), 

w — the true air speed, 

u« — the initial velocity relative to the air, 

*»> — the initial yaw relative to the morion of the 
air. 
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The initial velocity and yaw may be found by the 



formula* 

u,* = v # * -J- w* 2wv. sin Z cos A, (46) 

«.* a ( 1 - sin* Z cos* A ) w»/u.* (rad)* (47) 

In forward firing, A = 0 ; Z = 1600 mils, 
u* = v. + w, (46a) 

8, = 0. (47a) 

In rearward firing, A = 3200 mils, Z = 1600 
mils, 

u, = v. — w, (46b) 

8. = 0. (47b) 

The air density ratio H and the sound velocity 



ratio a/a. may be taken as constants. Approxi- 
mately, if Y is the altitude of the airplane above 
sea level, 

H s= e“«, h = 0.000,031,58 ft-*, (48) 

a/a. = e-«, k = 0.000,003,01 ft-*. (49) 

If s. is the standard stability factor, which cor- 
responds to the muzzle velocity, muzzle spin, and 
standard surface air density, the initial stability 
factor is 

*• = *»▼•*/**•* H, (50) 



If h is rhc yawing moment damping factor, 
% — the cross wind force damping factor, 

c'=(h + *)/2u, 
the loss in velocity due to yaw is 
Au _ 8A (s>— 1/2) A'u (u./a) 

a. (s.— 1) Cc' + c" (u*./a) 
Then the effective initial velocity is 



(51) 

(52) 



u. = u. — Au. 



(5?) 



Now we can find the distance measured along 
the line of departure to a point directly above the.* 
projectile, 

P = (C/H) [p(ua»/a) — p (ma./a)], (54) 

the drop of the projectile, 

Q = (Ca./Ha)* [q(ua./*> — q(u*./a)J 
- (Ca*/Ha) Pq, (u*./a), (55) 

and the time of flight is, 
t f = (Ca./Ha) [t(ua./a) — t (u^./a)]. (56) 

The horizontal and vertical components, if de- 
sired, may be found by the formulas 
x = Pcosfi* (57) 

y = P sin fi. — Q, (58) 

with the angle of departure determined by the 
formula 



% 

* (sinZ sin A)* + (sin Z cos A -f- w/v.) 1 

or 

cot* fi. = (tan Z sin A)* -f (tan Z cos A ■+■ w/v. cos Z)*. 



(59) 

(59a) 



Illustrative example is given on the following page. 
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(b) Example 
Data: 

Azimuth 
Zenith Angle 
Muzzle Velocity 
Air Speed 
Altitude 

Standard Stability Factor 
Damping Factors h a 

Ballistic Coefficient 



1200 mils 
800 mils 
3000 fps 
1000 ip* 
30,000 ft 
3.66 

1 0j6 sec- 1 
at 3000 fps 
0.235 on G|.i 



Problem: To find the time of flight and drop for 
a Siacd range (P) at 1200 ft 



Solution: (Use the Aircraft Table Based on 
Gs.i)« 

cos A = 0.38268 
sin Z = 0.70711 



sin Z cos A = 027060 



2wv, = 2x 1000x3000 = 6,000,000 



v,* = (3000)* 
w* = (1000)* 
2wv, sin Z cos A 
u#» 
u. 



= 9,000,000 
= 1000,000 
= 1,623,600 
= 11,623,600 (Eq. 46). 
= 3409.3 fps 



l*» = 0.72940 x 1,000,000/11,623,600 = 
0.06275 (Eq. 47) 



3. = 02505 rad (2552 mils), (1 rad = 
1018.6 mils). 

hY =0.000,031,58 x 30,000 = 0.9474 



H = 0.3877 (Eq. 48) 

kY = 0.000,003,01 x 30,000 = 0.0903 

a/a. = 0.9137 (Eq. 49) 

s. = 3.66 x 9,000,000/1 1,623,600 x 0.3877 = 
7.31 (Eq. 50) 

c' = 10.6/2 X 3000 = 0.00177 ft" 1 
ivu/a = 3409.3/0.91 37 = 373U fps 

A'u(<vu/a) = 1.4003 sec- 1 
(from Aircraft Table) 

c"(u^u/a) = 2288 x 10“* ft" 1 
(from Aircraft Table) 



u. = 3409.3 - 206.5 = 320Z8 fps (Eq. 53) 
Cd/H = 0235/0.3877 = 0^061 
Cc»a./Ha = 0.6061/0.9137 = 0.6633 
(Ccta./Ha)* = (0.6633)* =0.4400 
PH/Cti = 1200/0.6061 = 1979.9 ft 
u^./a = 3202.8/0.9137 = 3505.3 fps 
Find p(\gu/>) in Tahle, then p(tta«/a) by Eq. 
54. 



ua*/a p t q 

3174.3 15580.5 3.548 1657 £Ao6 

3505.3 13600.6 2.9$5 118.4 £>114 

1979.9 0.593 47.3 £0572 

t, = 0.6633 x 0.593 = 0.393 sec (Eq. 56) 

0.4400 x 47.3 = 20.8 

0.6633 x 1200x0.02114 = 168 
Q = 4.0ft (Eq. 55) 

(c) Interpolation. The derivatives listed in 
the Aircraft Tables should be used for inter- 
polating backward or forward from die nearest 
tabular value of the argument. For instance, in 
the last example, it is desired to find p for u = 

3505.3 and then to find u for p = 15580.5. First, 
the table gives 

at u= 3510, p = 13572.81, 

-^-= - 5.923 

du 

Hence, Au = —47, Ap = + 27.84, 

and, at u = 3505.3, p = 13600.65. 

Then we have p = 15580.5 
and the table gives 

p = 156062, - 6.044 at u = 3170. 

Hence, Ap = — 257 Au = -f 42, 

sod u= 31742. 



Au = 0.06275 x 0.9137 



6.81 x 1.4003 



6.31 x 0235 x 0.00177 + 0.000,022,88 

~ 0 05733 o m & TZ onw = 206 5 *• <*■ **> 
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b. Fin- stabilized Projectiles 

(1) Tabular Functions. Reference 6 contains 
a table of modified Siacci functions based on the 
urag coefficient for the 90-mm High Explosive 
Antitank Shell T108, which has rigid fins. These 
tables can be used for other similar projectiles, 
although they are not suitable for all fin-stabilized 
projectiles. The following functions and their 
derivatives are tabulated as functions of Mach 
number M : 



(60) 



(61) 



(62) 



(63) 



g/a, = 32.154/1 120.27 = 0.0287 see' 1 . 



Space 



Time 



S = 



T = 



f 

M 

2.7 

f- 

M 

2.7 



dM 

MKd(M) ’ 



dM 

M*Kd(M) * 



Inclination I 



Altitude 






dM 



M*K D (M) ’ 

* 

_ j' I(M)dM 



M 

2.7 



M 



MKd(M) ' 



The factor in the inclination function is 



y = x tan 8. — - ^ V^ *M A(M) - A(M.) ) 

tan * = tan - - **• - ?** — [I(M)-I(M.)J. 
The angle of departure for y = 0 may be found by 



sin 24, 



ao p, r A(M) - A(Mo) 
a *PY L S(M) - S(M 0 ) 




(2) Application 

(a) Formulas. The following formulas are 
applicable for computing trajectory data at super* 
elevations below 15°. 

Let 

m be the mass of the projectile in pounds, 
i — the form factor of the projectile, 

d — the maximum diameter of the projectile in 
feet, 

v — the velocity of the projectile in feet per 
second, 

v 0 — the muzzle velocity in feet per second, 
a — the velocity of sound in air in feet per second, 
p — the density of air in pounds per cubic foot, 
4 — the angle of inclination of the trajectory, 

4„ — the angle of departure. 

Also let 



a* = 1120.27 fps, 




p 0 = 0.07513 lb/ft 1 , 




M 0 = vo/a. 


(64) 


M == v cos $ /a cos I*. 


(65) 


T = p«id a /m. 


(66) 


Then the horizontal range in feet is 




* = (po/pr)lS(M) — S(M.)], 


(67) 



the time of flight in seconds is 
tf= (PosecVapT) [T(M)~T(M*)], (68) 



the altitude in feet is 



+ x 



aopo see 8 4» 
a * PT 



I(M*). 



the formula 



(69) 

(70) 

(71) 



If p = p 0 and a=ao, let S = S(M), S # =r 
S(M«), etc. Then these formulas can be ex- 
pressed more simply: 



* = (1/T)(S — So), 

tf = (sec. 4,/a,y) (T — To), 

y = x tan 4, — (A — Ao) +x 



(67a) 

(68a) 

SCC**. y 

a©T ** 
(69a) 



tan & = tan 4, 



sec 1 4. 



(I -I.). 



(:70a) 



and f or y = 0, 

< 71 *> 



14 




(b) Example 



Data: Gun, 90-mm ; Shell, HEAT, T106E20 



Caliber 

Maaa 

Muzzle Velocity 
Form Factor 
Air Density 
Sound Velocity 



029525 ft 
142 lb 
2400 fps 
1.00 on GnM 
0.07513 lb/ft* 
112027 fps 



Problem : To find the Te l o ci ty, time of flight, angle 
of departure, and angle of fall at a range 
of 3000 feet. 



Solution : From the Table of the Shod Functions 
baaed on Shell, 90-nun, HEAT, T408. 



Check : x tan 4, = 3000 (.00964) = 28.92 ft 
sec 3 4.(A - A*) 1.00010(0.0’.867) 

a.y* _ 0.000,23829 

= — 78.36 ft 



u= 

«-T 

= + 49.42 ft 



3000(1.00010) 

0.51667 



(0X0351) 



y = 28.92 + 49.42 - 78.36 = - 002 ft 
(Sq. 69a). 



7. DirFxaaNTiAi. Errxcrs 
Some of the effects on range and deflection 
caused by variations from standard conditions 
will be discussed and explained. Mon* detailed 
explanations are given in References 2 1 22, 23 
and 24. 



MS T A I 

123650 3X19196 1.427502 .02537 X) 1923 

2.14234* 120836 0209983 00670 00851 

Dfc 128360** 0217519 01867 X)1072 



* M. = 2400/112027 = 2.14234 (Eq. 64) 

T = 0.07513(1.00) (029525)*/142 = 
0.000,4612 ft" 1 (Eq. 66) 

a*T = 112027(0.000,4612) = 021667 sec" 1 

a»y* = 0.51667(0000,4612) = 

0.00023829 ft- 1 sec" 1 



**S- S. = T x = 0.000,4612(3000) = 
128360 (Eq. 67a) 



sin 24, = 



2 rOP1867 



0.51667 L 128360 
= 001928 (Eq. 71a) 



000851 



] 



24».= 19.64 mils, = 922 ««!«, cos 4, = 
0.99995, 



sec 4.= 1.00005, sec 1 4. = 100010 



tan 4' = 



100010 

0.51667 



(QXJ1G72) = 



000964 



— 0.02075 



tan 4= -0.01111 (Eq. 7th) 

4 = — 11 22 mils, cos 4 = 00994 
v = 112027 (173650) (0.99995) /(0.99994) 
= 1945.37 fps (Eq. 65) 

^S 0717519 ) = 1-3® «C (Eq. 68a) 



a. Effects on Range 

(1) Height of Target. Two methods may be 
used for determining the effect of the height of 
the target with respect to die gun : one by meant 
of the trajectory chart, and another by a correction 
to the angle of departure. 

(a) For large variations in height, die tra- 
jectory chart should be used. The trajectory that 
passes through the point determined by the g i ve n 
horizontal range and altitude can be found. The 
total range of this trajectory is the range where 
it crosses the axis. The diff erence between the 
given range and the total range is the eff ect of the 
height of target For greater accuracy, the ex- 
terior ballistic tables may be used in Ueu of the 
charts. 

(b) For small variations in height, an ap- 
proximate effect nay be found by adding the 
angle at site, 

« = tan- 1 (H/R'), (72) 

to the angle of departure f*. obtaining the cor- 
rected angle 

f = f. + «. (73) 

Here, H is the height of target, 

R< — the map range, 

— the angle of departure corresponding to 
the map range in the Exterior Ballistic 
Tables, 

f — the angle of departure correspaadmg to the 
co r rected range R. 

Then the effect on range due to height of target is 

A„R = R. — R. (74) 
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(2) Elevation. Since the Exterior Ballistic 
Tables give the range as a function of elevation for 

nitant values of muzzle velocity and ballistic 
coefficient, the effect of an increase in elevation 
may be found from the differences. At short 
ranges, the effect may be computed more ac- 
curately by using the Siacci tables, as explained 
above The vacuum formula (Eq. 6) may be use- 
ful in estimating the effect, espec ia lly for large 
ballistic coefficients. It should be noted that the 
effect is small in the vicinity of the maximum 
range and, of course, is negative at elevations 
above that corresponding to die maximum range. 

( 3) Muszle Velocity. Since the Exterior Bal- 
listic Tables give the range as a function of muzzle 
velocity for constant values of elevation and 
ballistic coefficient, the effect of an increase in 
muzzle velocity may be found fr o m the differences. 
At short ranges, the effect may be computed by 
using the Siacci tables. In vac uum , the range is 
proportional to the square of the velocity, as shown 
in Equation 6 ; in air, this is a good a p pro xim ation 
if the ballistic coefficient is large. 

( 4 ) Ballistic Coefficient. Since the Exterior 
Ballistic Tables give the range as a function of 
ballistic coefficient for constant values of muzzle 
velocity and elevation, the effect of an increase in 
log C may be found from the diffe rences. At 
short ranges, the effect may be computed by using 
the Siacci tables. It is customary to find the effect 
of a 1 % increase in C, which is equivalent to an 
increase of 0.00432 in log C. 

(5) Weight of Projectile. An increase m 
weight of projectile decreases the muzzle velocity 
and increases the ballistic coefficient. 

The effect on muzzle velocity may be ex pre ss ed 

A p v. = nv„(Ap)/p, (75) 

where v e is the muzzle v do a t y 
p is the projectile weight, 
n the logarithmic d iff ere nti al coefficient. 

The value of n may be oBtained from interior 
ballistic tables. It is ap pr ox imately — 0 3 for a 
rifled gun with multiperforated propellant grains, 
— 0.4 for a rifled gun with single-perfor a ted 
grains, —0.47 for a smooth bore mortar with 
flake propellant, and — 0.65 for a recoilless rifle 
with multiperforated grains. 

Since the ballistic co efficient is proportional to 
the projectile weight by definition. Equation (22), 



the effect on C may be expressed 

A,C = G(Ap)/p. (76) 

If A,R is the effect on range of 1 unit increase in 
muzzle velocity, 2«R the effect on range of 1% 
increase in ballistic coefficient, the effect on range 
of an increase in projectile weight is 

A*R = AyR(A^v) + AoR(100 A^Q/C (77) 

(6) Air Density. The retardation doe to air 
resistance is proportional to the air density and 
inversely proportional to the ballistic coefficient, 
as indicated in Eq. (25). The refo re, to a first 
approximation, the effect on range of an increase 
in air density is equal to the effect of the same 
percentage decrease in C 

(7) Air Temperature. An increase in the 
temperature of the atmosphere affects the r a ng e in 
two ways : it increases the density of the air and 
the velocity of sound. The temp era tur e of the 
air may affect the tem p erature of the propellant, 
which affects the muzzle velocity, but that is an 
interior ballistic problem. 

For a given pressure and humidity, the density 
is proportional to the absolute t empera ture. The 
effect of density on range has been treated above. 

For a given humidity, the velocity of sound is 
proportional to the square root of the absol u t e 
temperature. For a given projectile v eloc i ty, the 
Mach number is inversely proportional to the 
acoustic velocity. A change in M produces a 
change in the drag function. A new trajectory 
could be computed with the varied function ; how- 
ever, it is simpler to use an approximate formula 
deduced from a se mi -e mp irical formu la derived fay 
Dr. Gronwall :** 

A*R = 0.2929 [0.01R — 0.005v*Ml 

- ( I - 0.000,01 5y«)A©R] At, (78) 

where R is the range, 

AtR —the increase in range due to increase of 
temperature 

v, — the muzzle velocity in feet per second, 

A,R — die effect on R of 1 fps increase m v„ 

y « — die maximum ordinate in feet, 

AcR — die effect on R of 1% increase m C, 

At — the increase in t empera ture in de g r ees 
Faiirenhdt. 
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(8) Wind. A rear wind has three effects: it 
antes the coordinate system along with it rela- 
tive to the earth; it increases the angle of de- 
parture relative to the air; and it decreases the 
iwiwl» velocity relative to the air. The total 
gflect on range may be expressed as follows: 



A-R 



= w» [t + 



sin f AyR 
n Af 



f AvJ’ 
(79) 



where W* is the range co mp o n ent of the wind (a 
rear wind is p o siti ve; a bead wind is negative), 

T — the time of flight, 

f — the angle of departure, 

v* — the muzzle velocity, 

Af R — the effect on range doe to a change of 
Af radians in f. 

AvR — the effect on range due to a change Av. in 



b. Effects on Deflection 

(1) Wind. A crosswind has two effects: it 
moves the coordinate system with it and it changes 
the direction of motion relative to the air. The 
total e ff ect on linear deflection may be expressed 

A»Z = W» (T — R/v. cos f ), (80) 

where W B is the c ro ss w in d. 

The angular deflection due to wind, in radians, 

is 



Dw — W, (T/R — l/v # cos f). (81) 



(2) Drift. After the yaw doe to initial dis- 
turbances has damped out, a spinning shell has 
a yaw of repose due to the curva t ur e of the tra- 
jectory. This yaw is mainly to the right of the 
traj e ctory if the shell has right-handed spin; to 
the left if it has left-handed spin. It* may be de- 
termined at short ranges by firing alternately 
from barrels with right and left hand twists of 
rifling, or at long ranges by c or rect i ng the ob- 
served deflection for the effects of wind, cant of 
trunnions, and r otatio n of earth. 



The differential equation of motion ior the drift 
z is** 

- 2*gAvJLN . Gpx 

z mdnn*MN. x *C ' 

where g is the gravitational acceleration, 

A — the axial moment of inertia, 
v. — the muzzle velocity, 

L — the crosswind force (lift), 

N — the spin, 

N* — the initial spin, 
m — the mass of the projectile, 
d — the caliber, 

u — the velocity of the projectile relative to the 
air, 

M — the overturning moment about the center of 
gravity, 

x — the horizontal range, 

G — the drag function, 
p — the air density, 

p. — the standard air density (0.07513 lb/ft*), 

C — the ballistic coeffi ci e n t 

Dots denote derivatives, and the initial con- 
ditions are 
*• = 0, z. = 0. 

c. Ballistic Density, Temperature, and Wind 
In treating the effects of variations in the 
atmosphere, it has been assumed that die percent- 
age variations from standard density and tempera* 
tore and the wind are constant all along the tra- 
jectory. Actually, they vary with altitude, so 
that weighted mean values must be used. The 
weighting factors are determined by the propor- 
tion of the total effect that a variation occurring 
in each ahitnde zone would have. Instead of 
co m p u t in g a set of weightin g factors for each tra- 
jectory, a few r epres e nt a t ive curves are used. The 
sum of the products of the zonal variations by 
their weighting factors is called the ballistic 
density, temperature, or wind. The results are 
then treated as constants in computing the effects 
on range and deflecti on. 
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angle of departure, (elevation) As used in this 
pamphlet these terms mean quadrant angle of 
departure which is defined as the angle between 
the horizontal and the line along which the pro- 
jectile leaves the projector. One of the two 
initial conditions of a trajectory, die other con- 
dition being initial velocity. 

For a rocloet-fired projectile in which the pro- 
pellant burning continues after the rocket leaves 
the launcher, it is necessary to make assump- 
tions relative to the point at which the rocket 
b ecome s a missile in free flight and the velocity 
attained at that point 

ingle of inclination The angle between the hori- 
zontal and the tangent to the trajectory at any 
point At the origin of the trajectory this is the 
same as the angle of departure. 

apical (semi-apical) angle. In general the angle 
formed at the apex or tip of anything. As ap- 
plied to projectiles, the angle between the tan- 
gents to the curve outlining the contour of the 
projectile at its tip, or for semi-apical angle, 
the angle between the axis and one of the tan- 
gents. For a projectile having a conical tip, 
the cone apex angle. 

Dallistica. Branch of applied mechanics which deals 
with the motion and behavior characteristics of 
missiles, that is, projectiles, bombs, rockets, 
guided missiles, etc., and of accompanying phe- 
nomena. It can be conveniently divided into 
three branches : — interior ballistics, which deals 
with the motion of the projectile in the bore of 
the weapon ; exterior ballistics, which deals with 
the motion of the projectile while in flight ; and 
terminal ballistics, which is concerned with the 
effect and action of the projectile when it im- 
pacts or bursts. 

crosswind force. The component of the force of 
air resistance acting in a direction perpendicular 
to the direction of motion, and in the plane de- 
termined by the tangent to the trajectory formed 
by the center of gravity of a missile and the axis 
of the missile. (This plane is called the plane 
of yaw.) 

drag. The component of the force of air resistance 
acting in a direction opposite to the motion of 
the center of gravity of a missile. 



elevation. See angle of departure. 

exterior ballistics. See ballistics. 

horizontal range. The horizontal componou 
the distance from the beginning of the trajectory 
( poiat of origin) to any other point on the tra- 
jectory. 

initial velocity, (muzzle velocity) The projectile 
velocity at the moment that the projejr. t zeises 
to be acted upon by propelling forces. 

For a gun-fired projectile the initial velocity, 
expressed in feet or meters, per second, is called 
muzzle velocity or initial velocity. It is obtained 
by measuring the velocity over a distance for- 
ward of the gun, and correcting hack to the 
muzzle for the retardation in Bight. 

For a rocket-fired projectile, the velocity at- 
tained at the moment when rocket burnout oc- 
curs is the initial velocity. A slightly fictitious 
value is used, referred to an assumed point as 
origin of the trajectory. 

The initial velocity of a bomb dropped from 
an airplane is the speed of the airplane. 

Magnus force. The component of the force of air 
resistance acting in a direction perpendicular to 
the plane of yaw (see crosswind force). It is 
caused by the action between the boundary layer 
of air rotating with the projectile and the air 
stream. This force is small, compared with 
drag and crosswind force. 

muzzle velocity. The velocity with which the 
projectile leaves the muzzle of the gun. See 
initial velocity. 

slant range. The distance, measured along a 
straight line, from the beginning of the trajectory 
(point of origin) to a point on the trajectory. 

stability factor. A factor which indicates the rela- 
tive stability of a projectile under given condi- 
tions. It is dependent upon the force moments 
acting to align the projectile axis with the tan- 
gent to the trajectory, and the overturning force 
moments. For the projectile to be stable, the 
stability factor must be greater than unity. 

standard atmosphere. An average or representa- 
tive structure of the air used as a point of de- 
parture in computing Bring tables for projec- 
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tiles. The standard atmosphere adopted in 1956 
t the Armed Services and now known as the 
o. S. Standard Atmosphere, is that used by the 
ir.temational Civil Aviation Organization 
(ICAO). This Standard Atmosphere assumes 
a ground pressure of 760 m/m of mercury, a 
ground temperature of 15® C, air density at 
ground level of 0.076,475 fl>/ft\ and a tempera- 
ture gradient to the beginning of the strato- 
sphere expressed by the formula 

absolute temperature T(®K) =288.16-6J>H 

where H is height above sea level in k i l omete rs . 

In the stratosphere (above 11 kilometers) the 
temperature is assumed constant at 216l66® K_ 
For firing table purposes wind velocities of zero 
are assumed. 



supereittv.-tion. An added positive angle in anti- 
aircraft gunnery that compensates for the fall of 
the projectile during the time of flight due to the 
pull of gravity ; the angle the gun or projector 
must be elevated above the gun-target line. 

trajectory. The pith followed in spa ** by the 
center of gravity of a missile while in free flight 

yaw. The angle between the direction of motion 
of a projectile and the axis of the projectile, re- 
ferred to either as yaw, or more completely, 
angle of yaw. The angle of yaw increases with 
time of flight in an unstable projectile (see sta- 
bility factor ) and decreases to a constant value, 
called die yaw of repose , or the repost amgie of 
yaw, in a stable projectile. 
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Figure 1. Projectile Shapes 



21 





22 



MACH NO. 

Figure 2. Drag Coefficient vs Mach Number 
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Figure 6 (1 of 12). Trajectories for Projectile Type 2, 45° elevation, 1200 fps muznle velocity 
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Figure 6 (2 of 12). Trajectories for 
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Figure 6 (5 of 12). Trajectories for Projectile Type 2, 45* elevation, 4400 fps muzzle velocity 
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HORIZONTAL. RANG^-OOOO YDS.) 

Figure 6 (7 of 12). Trajectories for Projectile Type 2, 60° elevation, 1200 fps muzzle velocity 
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Figure 6 (8 of 12). Trajectories for Projectile Type 2, 60° elevation, 2000 ip e nr* 1 ' - '■ -* velo-it/ 
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Figure 
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HORIZONTAL (IOOO YDS) 

Figure 7 (4 of 12). Trajectoriea for Projectile Type 8, 45° elevation, 1^00 fps mu/wZie velocity 
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Figure 7 {5 of 12). Trajectories for Projectile Type 8, 45° elevation, 4400 fps* muzzle velocity 
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Figure 7 (Tof 12). Trajectories for Projectile Type 8, 60* elevation, 1200 fps mu Erie velocity 
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Figure 7 (9 of 12). Trajectories for Projectile Type 8 P 60* elevation, 2800 fps muzzle velocity 
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ENGINEERING DESIGN HANDBOOK SERIES 



The Engineering Design Handbook Series is intended to provide a compilation 
of principles and fundamental data to supplement experience in assisting engineers 
In the evolution of new designs which will meet tactical and technical needs while 
also embodying satisfactory producibility and maintainability. 

Listed below are the Handbooks which have been published or submitted for publi- 
cation. Handbooks with publication dates prior to 1 August 1962 were published as 
20-series Ordnance Corps pamphlets. AMC Circular 310-38, 19 July 1963, redesignated 
those publications as 706-series AMC pamphlets (e.g., ORDP 20-138 was redesignated 
AMCP 706-138). All new, reprinted, or revised handbooks are being published as 
706-series AMC pamphlets. 
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